**WALSH AND MCDONALD REPLICATION


****FIGURE 1****

*Figure 1a-1c


*Figure 1a: 2016 CCES
use "2016 CCES data.dta"

bysort troops: egen group_mean=wtmean(supportactnorm), weight(weight)
gen group_se=.
replace group_se=0.0160537 if troops==0
replace group_se=0.017246 if troops==1
gen conf_high = group_mean + 1.96*group_se
gen conf_low = group_mean - 1.96*group_se
gen Z_jitter = troops + runiform(-0.3, 0.3)
gen Y_jitter = supportactnorm + runiform(-0.05, 0.05)

#delimit ;
twoway 
	scatter Y_jitter Z_jitter, mcolor(%10) msize(1) || 
	scatter group_mean troops, msize(1) mcolor(black) || 
	rcap conf_high conf_low  troops, color(black);

clear all	


*Figure 1b: 2018 CCES
use "2018 CCES data.dta"

bysort troops: egen group_mean=wtmean(supportactnorm), weight(teamweight)
gen group_se=.
replace group_se=0.0168858 if troops==0
replace group_se=0.016158 if troops==1
gen conf_high = group_mean + 1.96*group_se
gen conf_low = group_mean - 1.96*group_se
gen Z_jitter = troops + runiform(-0.3, 0.3)
gen Y_jitter = supportactnorm + runiform(-0.05, 0.05)

#delimit ;
twoway 
	scatter Y_jitter Z_jitter, mcolor(%10) msize(1) || 
	scatter group_mean troops, msize(1) mcolor(black) || 
	rcap conf_high conf_low  troops, color(black);

clear all


*Figure 1c: 2019 MTurk
use "2019 MTurk.dta"

bysort troops: egen group_mean=mean(supportactnorm)
gen group_se=.
replace group_se=0.0130052 if troops==0
replace group_se=0.0139343 if troops==1
gen conf_high = group_mean + 1.96*group_se
gen conf_low = group_mean - 1.96*group_se
gen Z_jitter = troops + runiform(-0.3, 0.3)
gen Y_jitter = supportactnorm + runiform(-0.05, 0.05)

#delimit ;
twoway 
	scatter Y_jitter Z_jitter, mcolor(%10) msize(1) || 
	scatter group_mean troops, msize(1) mcolor(black) || 
	rcap conf_high conf_low  troops, color(black);

clear all

*Meta-Analysis to go with Fig. 1
use "C:\Users\Jared\Dropbox\Drones\metaanalysis.dta"
metan effect lowerbound upperbound, fixed

clear all


****FIGURE 2****
use "2019 MTurk.dta"

bysort troops: egen group_mean=mean(milcasnorm)
gen group_se=.
replace group_se=0.014276 if troops==0
replace group_se=0.0126973 if troops==1
gen conf_high = group_mean + 1.96*group_se
gen conf_low = group_mean - 1.96*group_se
gen Z_jitter = troops + runiform(-0.3, 0.3)
gen Y_jitter = milcasnorm + runiform(-0.05, 0.05)

#delimit ;
twoway 
	scatter Y_jitter Z_jitter, mcolor(%10) msize(1) || 
	scatter group_mean troops, msize(1) mcolor(black) || 
	rcap conf_high conf_low  troops, color(black);



****ONLINE APPENDIX****


*Table A1: Effect of Troops Treatment on Attack Support (Weighted where weights were available)

*2016 CCES (Weighted)
use "2016 CCES data.dta"
reg supportactnorm troops [aw=weight], vce(hc2)
*Minimum Detectable Effect
power twomeans 0.767116, sd(.257) power(0.8) n(805)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
reg supportactnorm troops [aw=teamweight], vce(hc2)
*Minimum Detectable Effect
power twomeans 0.6087743, sd(.295) power(0.8) n(879)

*2019 MTurk (Unweighted)
use "2019 MTurk.dta"
reg supportactnorm troops, vce(hc2)
*Minimum Detectable Effect
power twomeans 0.6052539, sd(.289) power(0.8) n(922)



*Table A2: Effect of Troops Treatment on Attack Support (Unweighted)

*2016 CCES (Weighted)
use "2016 CCES data.dta"
reg supportactnorm troops, vce(hc2)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
reg supportactnorm troops, vce(hc2)

*2019 MTurk (Unweighted)
use "2019 MTurk.dta"
reg supportactnorm troops, vce(hc2)



*Table A3: Effect of Troops Treatment on Attack Support (Weighted, with covariates)

*2016 CCES (Weighted)
use "2016 CCES data.dta"
reg supportactnorm troops white age male faminc ideo5 [aw=weight], vce(hc2)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
reg supportactnorm troops white age male faminc ideo5 [aw=teamweight], vce(hc2)

*2019 MTurk (Unweighted)
use "2019 MTurk.dta"
reg supportactnorm troops white age male faminc ideo5, vce(hc2)



*Table A4: Effect of Troops Treatment on 1) Support for Attack, and 2) Expectation of Military Casualties

use "2019 MTurk.dta"
*Support for Attack
reg supportactnorm troops, vce(hc2)
reg supportactnorm drones troops, vce(hc2) nocons
*Expected Military Casualties
reg milcasnorm troops, vce(hc2)
reg milcasnorm drones troops, vce(hc2) nocons


*Table A5: Distribution of Responses on Expected Military Casualties by Condition

use "2019 MTurk.dta"
tab milcasnorm troops, col



*Table A6: Distribution of Responses on Expected Civilian Casualties by Condition

use "2019 MTurk.dta"
tab civcasnorm troops, col



*Table A7: Experimental Results by Veteran Status (Weighted, Robust Standard Errors)

*2016 CCES (Weighted)
use "2016 CCES data.dta"
bysort milstat_3: reg supportactnorm troops [aw=weight], vce(hc2)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
bysort milstat_3: reg supportactnorm troops [aw=teamweight], vce(hc2)



*Table A8: Experimental Results by Partisanship (Weighted, Robust Standard Errors)

*2016 CCES (Weighted)
use "2016 CCES data.dta"
bysort pidwleaners: reg supportactnorm troops [aw=weight] if pidwleaners!=2, vce(hc2)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
bysort pidwleaners: reg supportactnorm troops [aw=teamweight] if pidwleaners!=2, vce(hc2)



*Table A9: Randomization Checks - Predicting Assignment to Condition Using Sociopolitical Characteristics

*2016 CCES (Weighted)
use "2016 CCES data.dta"
probit troops white age male faminc ideo5 [pw=weight]

*2018 CCES (Weighted)
use "2018 CCES data.dta"
probit troops white age male faminc ideo5 [pw=teamweight]

*2019 MTurk (Weighted)
use "2019 MTurk.dta"
probit troops white age male faminc ideo5



*Table A10: Descriptive Statistics across All Three Studies – Standardized Differences in Means

*2016 CCES (Weighted)
use "2016 CCES data.dta"
sum supportactnorm if drones==1 [aw=weight]
sum supportactnorm if troops==1 [aw=weight]
*Calculating Cohen's D
esizei 413  .767116 .2573997    392 .7599028 .2571285
*Equivalence Test
tostt supportactnorm, by(troops) eqvtype(delta) eqvlevel(0.08)

*2018 CCES (Weighted)
use "2018 CCES data.dta"
sum supportactnorm if drones==1 [aw=teamweight]
sum supportactnorm if troops==1 [aw=teamweight]
*Calculating Cohen's D
esizei 439  .6087743 .2995329    440 .6265754 .2929721
*Equivalence Test
tostt supportactnorm, by(troops) eqvtype(delta) eqvlevel(0.08)

*2019 MTurk (Unweighted)
use "2019 MTurk.dta"
sum supportactnorm if drones==1
sum supportactnorm if troops==1
*Calculating Cohen's D
esizei 445  .6052539 .2939566    477 .5782788 .2840268
*Equivalence Test
tostt supportactnorm, by(troops) eqvtype(delta) eqvlevel(0.08)


